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ABSTRACT 


This  report  deals  with  the  design  and  performance  evaluation  of  a  distributed  de¬ 
tection  system  for  a  dispersing  biological  warfare  (BW)  cloud  embedded  in  the  natural 
aerosol  component  of  the  background  air.  The  distributed  system  employs  a  number  of 
physically  separated  BW  agent  point  detectors  (sentries)  located  within  some  target  region 
and  a  data  fusion  center  that  provides  the  final  decision  as  to  the  presence  or  absence  of 
the  bio-target  by  combining  the  individual  localized  decisions  from  the  various  point  de¬ 
tectors  using  a  prespecified  combining  strategy.  In  this  system,  each  detector  implements 
a  generalized  likelihood  ratio  test  on  its  own  localized  observations  to  test  for  the  presence 
or  ab.sence  of  a  bio-target.  These  localized  detection  decisions  are  then  transmitted  to  a 
data  fusion  center  where  they  are  logically  combined  to  yield  a  global  detection  decision 
for  the  distributed  system.  The  optimization  of  the  global  detection  performance  of  the 
distributed  system  is  derived  by  application  of  the  Lagrange  multipliers  method,  whereby 
the  global  probability  of  detection  for  the  system  is  maximized  subject  to  the  constraint 
that  the  global  probability  of  false  alann  is  maintained  at  a  prespecified  level  (viz.,  at  a 
constant  and  tolerable  false-alarm  rate).  A  number  of  combining  strategies  are  investigated 
in  order  to  determine  some  overall  .system  optimality  for  detection.  It  is  found  that  the 
optimal  detection  strategy  for  each  individual  BW  agent  detector  in  the  distributed  system 
depends  on  the  strategy  of  all  the  other  detectors,  as  well  as  on  the  structure  of  the  data 
fusion  center.  Explicit  numerical  results,  in  the  form  of  the  probability  of  detection  versus 
signal-to-noise  ratio  for  several  preassigned  false-alarm  probabilities  and  sample  sizes  [which 
determine  the  maximum  allowable  mean  time  to  detection],  are  presented  for  the  case  of 
a  distributed  system  consisting  of  three  BW  agent  detectors.  The  results  are  compared 
to  the  detection  performance  that  is  achievable  using  only  a  single  detector.  Finally,  the 
problem  of  the  correct  placement  of  the  individual  BW  agent  detectors  in  the  distributed 
system  in  order  to  achieve  a  specified  detection  performance  is  addressed.  In  particular,  a 
method  for  the  determination  of  the  dosage  probabilities  at  fixed  points  in  the  target  region 
is  developed  and  used  to  place  BW  agent  detectors  in  a  distributed  system. 
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EXECUTIVE  SUMMARY 


Title:  Methodology  for  the  Design  and  Optimal  Placement  of  Point  Detectors  in  a  Dis¬ 
tributed  Detection  System  for  Remote  Defence  Against  Biological  Warfare  Agent  Releases 

Author:  Eugene  Yee 

Introduction:  There  are,  quite  rightly,  growing  concerns  world-wide  about  the  dangers, 
both  actual  and  potential,  of  the  use  of  biological  warfare  (BW)  agents.  In  particular,  the 
May  1997  Report  of  the  Quadrennial  Defence  Review  (Department  of  Defence,  USA) 
concluded  that  the  likelihood  of  encountering  BW  agents  in  a  specific  military  operation 
increases  with  each  year  and  will  probably  be  a  likely  condition  of  future  warfare.  The  low 
cost,  low  technology,  and  relative  ease  of  producing  BW  agents,  and  our  national  policy  to 
maintain  appropriate  defences  against  such  threats  imply  that  the  BW  agent  threat  must 
be  taken  seriously.  Indeed,  the  use  of  BW  agents  against  Canadian  forces  is  probably  more 
likely  than  it  was  in  the  past  given  the  proliferation  of  BW  agent  capabilities  among  various 
countries. 

The  lack  of  direct  BW  agent  experience  implies  that  quantified  military  hazard 
assessment  of  BW  agents  must  be  approached  primarily  as  a  description  of  possibilities, 
largely  based  on  considerations  of  past  research  and  development  activities  in  this  area 
and  supported  by  mathematical  modelling  of  the  possible  risks  to  an  individual  or  a  troop 
of  being  exposed  to  concentrations  of  a  dispersing  BW  agent  cloud  released  from  some 
form  of  delivery/munition  system.  Depending  on  the  release  mode,  BW  agent  sources  can 
include  transient  line  and  point  sources  produced  by  spraying  from  a  low-flying  aircraft 
or  some  stationary  platform,  as  well  as  multiple  instantaneous  point  sources  produced  by 
sub-munitions  and  bomblets.  Furthermore,  BW  agents  are  extremely  potent,  and  attacks 
that  are  initiated  by  disseminating  relatively  small  quantities  of  the  agent  either  overtly  or 
covertly  from  a  fixed  or  mobile  platform  .such  as  a  truck,  ship,  or  fixed-wing  aircraft  can 
produce  a  significant  hazard  range.  In  consequence,  in  order  to  develop  an  effective  defence 
against  BW  agent  attack,  it  is  necessary  to  understand  the  course  and  consequences  of 
such  an  attack  ,  so  that  the  design  and  evaluation  of  detection  and  warning  systems  can  be 
facilitated.  The  prediction  of  what  will  happen  to  a  BW  agent  after  it  has  been  disseminated 
from  a  munition  system  and  the  evaluation  of  expected  concentration  levels  available  for 
BW  agent  detection  systems  located  downwind  of  the  release  position  can  only  be  obtained 
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using  computer  simulation  models. 

When  a  BW  agent  is  released  into  the  atmo.sphere,  the  immediate  fate  of  the  agent 
largely  depends  on  the  local  meteorological  conditions.  In  almost  all  cases,  drift  and  tur¬ 
bulent  diffusion  will  spread  the  BW  agent  cloud  over  a  large  area.  Depending  on  the  at¬ 
mospheric  turbulence,  the  BW  agent  cloud  will  become  more  or  less  dispersed.  However,  it 
must  be  recognized  that  local  micrometeorology  can  easily  dominate  the  drift  and  diffusion 
process.  The  turbulence  is  dependent  on  the  ground  topography,  the  temperature  in  the  air 
mass,  the  wind  speed,  etc.  and  the  actual  distribution  and  concentration  of  a  BW  agent  in  a 
dispersing  cloud  can  vary  greatly  from  place  to  place.  Any  BW  agent  detection  and  warning 
system  must  be  designed  to  recognize  and  account  for  these  possibilities.  In  view  of  this,  the 
theory  of  distributed  detection  of  BW  agents  is  developed  in  this  report.  Furthermore,  it 
is  shown  how  a  mathematical  model  for  simulating  the  turbulent  dispersion  of  BW  agents 
in  the  atmosphere  can  be  applied  to  formulate  a  rational  and  practical  methodology  for 
the  optimal  placement  of  individual  BW  agent  detectors  in  a  remote  distributed  detection 
system. 


Results:  This  report  deals  with  the  design  and  performance  evaluation  of  a  distributed 
detection  system  for  a  dispersing  biological  warfare  (BW)  cloud  embedded  in  the  natural 
aerosol  component  of  the  background  air.  The  distributed  system  employs  a  number  of 
physically  separated  BW  agent  point  detectors  (sentries)  located  within  some  target  region 
and  a  data  fusion  center  that  provides  the  final  decision  as  to  the  presence  or  absence 
of  the  bio-target  by  combining  the  individual  localized  decisions  from  the  various  point 
detectors  using  a  prespecified  combining  strategy.  In  this  system,  each  detector  implements 
a  generalized  likelihood  ratio  test  on  its  own  localized  observations  to  test  for  the  presence 
or  absence  of  a  bio-target.  These  localized  detection  decisions  are  then  transmitted  to  a 
data  fusion  center  where  they  are  logically  combined  to  yield  a  global  detection  decision 
for  the  distributed  system.  The  optimization  of  the  global  detection  performance  of  the 
distributed  system  is  derived  by  application  of  the  Lagrange  multipliers  method,  whereby 
the  global  probability  of  detection  for  the  system  is  maximized  subject  to  the  constraint 
that  the  global  probability  of  false  alarm  is  maintained  at  a  prespecified  level  (viz.,  at  a 
constant  and  tolerable  false-alarm  rate).  A  number  of  combining  strategies  are  investigated 
in  order  to  determine  some  overall  system  optimality  for  detection.  It  is  found  that  the 
optimal  detection  strategy  for  each  individual  BW  agent  detector  in  the  distributed  system 
depends  on  the  strategy  of  all  the  other  detectors,  as  well  as  on  the  structure  of  the  data 
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fusion  center.  Explicit  numerical  results,  in  the  form  of  the  probability  of  detection  versus 
signal- to-noise  ratio  for  several  preas.signed  false-alarm  probabilities  and  sample  sizes  [which 
determine  the  maximum  allowable  mean  time  to  detection],  are  presented  for  the  case  of 
a  distributed  system  consisting  of  three  BW  agent  detectors.  The  results  are  compared 
to  the  detection  performance  that  is  achievable  using  only  a  single  detector.  Finally,  the 
problem  of  the  correct  placement  of  the  individual  BW  agent  detectors  in  the  distributed 
system  in  order  to  achieve  a  specified  detection  performance  is  addressed.  In  particular,  a 
method  for  the  determination  of  the  dosage  probabilities  at  fixed  points  in  the  target  region 
is  developed  and  used  to  place  BW  agent  detectors  in  a  distributed  system. 


Significance  of  Results:  The  design  and  optimal  placement  of  BW  agent  detectors  in 
a  distributed  system  for  the  space-time  detection  of  weak  target  bio-aerosol  signals  from 
a  dispersing  BW  agent  cloud  embedded  in  a  general  aerosol  background  clutter  have  been 
investigated.  The  results  of  this  report  provide  a  general  framework  for  obtaining  potential 
improvements  of  distributed  systems  for  the  remote  detection  of  BW  agent  releases  by  using 
appropriate  fusion  rules  or  combining  strategies  at  the  data  fusion  center.  These  combining 
strategies  can  be  used  to  improve  the  probability  of  detection  while  maintaining  a  small 
(and,  hence,  tolerable)  false-alarm  rate,  and  coupled  with  a  method  for  the  correct  placement 
of  the  individual  BW  agent  detectors  of  the  distributed  system  in  the  target  region,  it  is 
now  possible  to  maximize  detection  performance  for  a  given  detection  environment  and  BW 
agent  scenario.  An  innovation  of  the  current  methodology  is  the  use  of  a  state-of-the-art 
atmospheric  dispersion  model  to  determine  the  area  covered  by  a  certain  critical  dosage 
with  a  given  probability,  and  to  use  this  information  for  the  optimal  placement  of  BW 
agent  detectors  in  a  distributed  system.  In  other  words,  a  transport  and  diffusion  model  for 
atmospheric  dispemion  is  used  to  quantify  the  detection  enviroment,  and  this  information 
is  then  used  to  obtain  explicit  solutions  for  optimum  thresholds  and  strategies. 

Of  particular  significance  is  the  fact  that,  for  a  given  set  of  operating  conditions 
(e.g.,  signal-to-noise  ratio,  sample  size,  false-alarm  probability,  etc.),  the  distributed  system 
always  exhibits  better  detection  performance  than  that  afforded  by  a  single  detector  system. 
Hence,  in  the  optimal  remote  detection  of  BW  agents,  it  is  preferable  to  have  many  low-cost 
BW  agent  detectors  configured  to  act  as  a  distributed  system,  rather  than  a  few  high-cost 
BW  agent  detectors  each  of  which  is  configured  to  act  in  a  stand-alone  mode.  Furthermore, 
distributed  detection  systems  ai-e  expected  to  be  more  reliable  and  immune  to  background 
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clutter  than  stand-alone  systems.  While  these  properties  are  due  to  inherent  redundancy  in 
a  distributed  system,  the  enhanced  performance  of  distributed  systems  essentially  derives 
from  the  coverage  and  spatial-temporal  diversity  afforded  by  the  system  configuration.  Fi¬ 
nally,  within  the  general  framework  of  distributed  detection  systems,  our  analysis  establishes 
that  the  optimal  decision  strategy  of  any  BW  agent  detection  system  must  take  into  ac¬ 
count  a  number  of  related  factors:  namely,  the  correct  placement  of  the  individual  detectors 
within  the  target  region  and  the  recognition  that  even  with  the  correct  placement  of  detec¬ 
tors,  the  optimum  decision  strategy  of  any  detector  must  take  into  account  the  operating 
environments  and  performances  of  all  other  detectors  in  the  system,  as  well  as  the  structure 
of  the  fusion  center. 

Future  Work:  Future  work  on  the  distributed  detection  of  BW  agents  is  possible  in  sev¬ 
eral  directions.  Firstly,  it  is  recommended  that  the  proposed  methodology  for  the  design 
and  placement  of  BW  agent  detectors  in  a  distributed  system  be  integi'ated  with  a  practi¬ 
cal  probabilistic  model  for  turbulent  diffusion  such  as  the  Second-order  Closure  Integrated 
PUFF  (SCIPUFF)  model  developed  by  Aeronautics  Research  Associates  Princeton  (ARAP) 
under  the  sponsorship  of  the  U.S.  Defence  Special  Weapons  Agency  (DSWA).  Secondly,  it 
would  be  meaningful  to  search  for  sub-optimum,  and  possibly  robust  designs  involving  far 
less  complex  computations  than  tho.se  proposed  for  the  current  methodology.  Finally,  it 
would  be  interesting  to  look  at  different  topologies  of  distributed  BW  agent  detection  sys¬ 
tems,  and  investigate  whether  certain  topologies  are  more  appropriate  for  use  in  defending 
against  certain  BW  agent  release  scenarios. 
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I.  INTRODUCTION 

There  axe,  quite  rightly,  growing  conceras  world-wide  about  the  dangers,  both  actual 
and  potential,  of  the  use  of  biological  warfare  (BW)  agents.  In  particular,  the  likelihood  of 
encountering  BW  agents  in  a  specific  military  operation  increases  with  each  year  [1,2].  The 
low  cost,  low  technology,  and  relative  ease  of  producing  BW  agents,  and  our  national  policy 
to  maintain  appropriate  defences  against  such  threats  imply  that  the  BW  threat  must  be 
taken  seriously.  Indeed,  the  use  of  BW  agents  against  Canadian  forces  is  probably  more 
likely  than  it  was  in  the  past  given  the  proliferation  of  BW  agent  capabilities  among  various 
countries. 

The  lack  of  direct  BW  agent  experience  implies  that  quantified  military  hazard 
assessment  of  BW  agents  must  be  approached  primarily  as  a  description  of  possibilities, 
largely  based  on  considerations  of  past  research  and  development  activities  in  this  area  and 
supported  by  mathematical  modelling  of  the  possible  risks  to  an  individual  or  a  troop  of 
being  exposed  to  concentrations  of  a  dispersing  BW  agent  cloud  released  from  some  form 
of  delivery/munition  system.  Depending  on  the  release  mode,  the  BW  agent  sources  can 
include  transient  line  and  point  sources  (e.g.,  spraying  from  a  low-fiying  aircraft  or  some 
stationary  platform)  as  well  as  multiple  instantaneous  point  sources  (e.g,  sub-munitions  and 
bomblets).  Furthermore,  BW  agents  are  extremely  potent,  and  attacks  that  are  initiated  by 
disseminating  relatively  small  quantities  of  the  agent  either  overtly  or  covertly  from  a  fixed 
or  mobile  platform  (e.g.,  truck,  ship,  fixed-wing  aircraft,  etc.)  can  produce  a  significant 
hazard  range  (i.e.,  the  maximum  downwind  distance  that  might  be  reached  by  a  militarily 
significant  concentration  of  agent  over  a  given  exposure  time  to  produce  an  infective  or 
effective  dose  under  particular  release  conditions).  In  consequence,  in  order  to  develop  an 
effective  defence  against  a  BW  agent  attack,  it  is  necessary  to  understand  the  course  and 
consequences  of  such  an  attack  (e.g.,  transport  and  degree  of  contamination)  so  that  the 
design  and  evaluation  of  detection  and  warning  systems  can  be  facilitated.  The  prediction  of 
what  will  happen  to  a  BW  agent  after  it  has  been  disseminated  from  a  munition  system  and 
the  evaluation  of  expected  concentration  levels  available  for  BW  agent  detection  systems 
located  downwind  of  the  release  position  can  only  be  obtained  using  computer  simulation 
models. 


When  a  BW  agent  is  released  into  the  atmosphere,  the  immediate  fate  of  the  agent 
largely  depends  on  the  local  meteorological  conditions.  In  almost  all  cases,  drift  and  tur- 
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bulent  diffusion  will  spread  the  BW  agent  cloud  over  a  large  area.  Depending  on  the  at¬ 
mospheric  turbulence,  the  BW  agent  cloud  will  become  more  or  less  dispersed.  However,  it 
must  be  recognized  that  local  micrometeorology  can  easily  dominate  the  drift  and  diffusion 
process.  The  turbulence  is  dependent  on  the  ground  topography,  the  temperature  in  the  air 
mass,  the  wind  speed,  etc.  and  the  actual  distribution  and  concentration  of  a  BW  agent  in  a 
dispersing  cloud  can  vary  greatly  from  place  to  place.  Any  BW  agent  detection  and  warning 
system  must  be  designed  to  recognize  and  account  for  these  possibilities.  In  view  of  this,  the 
theory  of  distributed  detection  of  BW  agents  is  developed  in  this  report.  Furthermore,  it 
is  shown  how  a  mathematical  model  for  simulating  the  turbulent  dispersion  of  BW  agents 
in  the  atmosphere  can  be  applied  to  formulate  a  rational  and  practical  methodology  for 
the  optimal  placement  of  individual  BW  agent  detectors  in  a  remote  distributed  detection 
system. 

The  consideration  of  a  distributed  system  as  opposed  to  the  single  detector  scheme 
for  BW  agent  detection  is  motivated  primarily  by  the  potentially  improved  detection  perfor¬ 
mance  at  a  preassigned  false-alarm  probability  that  can  be  achieved  using  multiple  detectors. 
Indeed,  the  performance  of  a  single  sensor  detection  system  is  rather  limited  because  the 
perfoi-mance  of  the  sensor  is  based  on  its  local  operating  conditions.  Consequently,  the  use 
of  a  distributed  sensor  detection  scheme  can  make  the  detection  performance  robust  in  ad¬ 
dition  to  providing  an  improvement  in  the  detection  probability  for  a  fixed  false-alarm  rate. 
However,  the  potential  improvement  in  detection  performance  achievable  in  a  distributed 
system  requires  a  careful  examination  of  two  factors;  namely,  (1)  the  design  of  the  data 
processing  scheme  used  to  proce.ss  the  observations  received  by  the  various  detectors;  and, 
(2)  the  development  of  a  methodology  for  the  optimal  placement  of  the  individual  detectors 
of  the  distributed  system  within  the  target  region  in  order  to  account  for  the  influences  of  a 
number  of  meteorological  parameters  (wind  speed,  turbulence,  topography,  etc.)  that  deter¬ 
mine  the  precise  distribution  and  concentration  of  a  dispersing  BW  agent  cloud.  Neglecting 
either  one  of  these  two  factors  can  lead  to  a  serious  degradation  in  the  overall  performance 
of  the  BW  agent  detection  system. 

This  study  focusses  on  the  formulation,  design,  and  analysis  of  a  distributed  and 
decentralized  detection  scheme  for  a  BW  agent  that  has  been  transported  from  some  re¬ 
mote  upwind  source  location  over  the  receptor  point  by  the  movements  of  the  wind.  In  this 
approach,  a  number  of  peripheral  BW  agent  detectors  (e.g.,  sentries)  are  deployed  over  a 
potential  target  area.  The  individual  detectors  are  assumed  to  process  their  local  observa¬ 
tions  through  some  optimal  threshold  detection  algorithm  (e.g.,  a  likelihood  ratio  detector) 
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to  produce  their  own  local  detection  decisions.  The  detection  decisions  from  the  individual 
detectors  are  then  transmitted  to  the  data  fusion  center  which  then  logically  combines  these 
local  decisions  to  decide  on  the  presence  or  absence  of  the  BW  agent.  Within  this  frame¬ 
work,  we  study  the  optimal  design  for  combining  decisions  emanating  from  peripherally 
located  detectors.  The  principal  objective  is  to  formulate  a  globally  optimum  detection  al¬ 
gorithm  for  the  distributed  system  as  a  whole,  in  terms  of  a  number  of  relevant  parameters; 
namely,  the  probability  of  detection  for  a  preassigned  false- alarm  probability,  the  signal- to- 
noise  ratio  presented  by  the  operating  environment,  the  minimum  mean  time  to  detection, 
and  the  location  or  placement  of  the  individual  detectors  or  sentries  within  the  designated 
target  area.  The  principal  task  in  this  report  is  to  combine  all  these  relevant  parameters 
to  obtain  the  “macro- algorithms”,  so  to  speak,  for  the  entire  distributed  BW  agent  detec¬ 
tion  system,  and  to  analyze  their  expected  performance  measures  (i.e.,  here  probabilities  of 
correct  decisions)  in  such  a  way  as  to  incorporate  the  controlling  physical  factors  which  the 
detection  environment  and  release  scenario  impose  (and  which  is  characterized  here  by  the 
application  of  mathematical  models  for  the  prediction  and  assessment  of  BW  agent  hazards 
in  the  atmosphere). 

II.  DETECTION  PERFORMANCE  OF  A  SINGLE  DETECTOR 


The  detectors  (sentries)  in  the  distributed  system  are  assumed  to  act  independently 
of  each  other.  In  order  to  characterize  the  detection  performance  of  the  distributed  system, 
it  is  necessary  to  first  model  the  detection  performance  for  a  single  detector.  To  proceed 
further,  it  is  necessary  to  either  measure  the  detection  performance  of  the  detector  as  a 
function  of  the  signal-to-noise  ratio  of  the  operating  environment,  or  to  develop  a  theoretical 
model  that  can  be  used  to  predict  the  detection  performance.  In  this  section,  we  do  not 
present  an  analysis  of  any  specific  BW  agent  detector.  Rather,  we  will  characterize  the 
detection  performance  of  a  generic  BW  agent  detector  in  order  to  use  it  as  an  example  of 
the  proposed  methodology  for  distributed  system  design  and  detector  placement. 

To  this  purpose,  consider  a  generic  BW  agent  detector  which  samples  for  a  total 
duration  T  to  produce  count  data  as  a  set  of  numbers  {F{j),j  =  1,2, . . . ,  J}.  Normally, 
the  count  data  is  related  to  the  time-averaged  concentration  measured  over  the  sampling 
interval  (e.g.,  the  count  in  each  sample  interval  is  proportional  to  the  number  of  particles, 
AT,  in  a  given  volume,  Rt,  of  air  where  R  is  the  known  rate  at  which  air  is  drawn  through  the 
detector  and  r  is  the  sampling  interval,  with  the  constant  of  proportionality  here  dependent 
on  a  correction  for  the  efficiency  of  collection,  retention,  or  absorption).  The  count  data 
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corresponding  to  the  target  signal  can  be  viewed  as  being  proportional  to  the  dosage  V 
accumulated  by  the  detector  during  the  time  r,  viz. 

Fiti)  =  Fij)  oc  r)=  x{t')  dt',  (1) 

Jtj 

where  xit)  is  the  instantaneous  concentration  (number  of  particles  per  unit  volume)  at 
time  t.  Furthermore,  the  count  data  sampled  by  the  detector  are  assumed  to  depend  lin¬ 
early  on  the  characteristic  signal  that  is  used  to  detect  the  bio-aerosol  (e.g.,  fluorescence 
energy,  etc.).  As  an  example,  in  laser-induced  fluorescence  techniques  for  BW  agent  detec¬ 
tion,  the  count  data  would  be  related  directly  to  bio-fluorescence  energy  (e.g.,  number  of 
fluorescing  photons)  in  a  characteristic  emission  band  [3] .  The  sample  counts  are  random 
variables.  It  is  assumed  that  the  count  F{j)  obtained  during  the  jth  sampling  interval  (re¬ 
ception  or  snapshot)  is  a  Poisson  random  variable.  Furthermore,  it  is  assumed  that  counts 
obtained  during  disjoint  sampling  intervals  (i.e.,  exposures)  are  mutually  independent  ran¬ 
dom  variables. 


The  count  data  are  assumed  to  consist  of  the  sum  of  a  target  bio-aerosol  signal 
(e.g.,  BW  agent)  and  an  ambient  interfering  background  clutter  aerosol  signal  (e.g.,  non- 
microbiological  aerosols,  pollen,  dust,  etc.)  that  is  unavoidably  received  from  the  natural 
aerosol  component  present  in  the  atmosphere.  In  this  case,  the  count  datum  F{j)  is  a 
PoLsson  random  variable  with  parameter  =  S(j)  -1-  For  simplicity,  we  take  the 

target  signal  S{j)  to  be  equal  to  the  partial  dosage  from  the  BW  agent  cloud  accumulated 
during  the  sampling  interval  r  extending  from  tj  to  tj  4-  t  (i.e.,  we  ignore  the  collection 
efficiency  of  the  sampler) ,  so 


S{j)  =  Vij)  =V{tf,r)  = 


dt\ 


where  Xsit)  is  the  instantaneous  concentration  contributed  by  the  BW  agent  cloud  (if 
present).  We  assume  that  the  target  signal  contribution,  S{j),  to  the  Poisson  parameter 
A(^)(j)  is  known  a  priori.  In  addition,  it  is  assumed  that  the  background  clutter  contri¬ 
bution,  A^°\  to  the  Poisson  parameter  X^^\j)  does  not  depend  on  j,  viz.,  the  background 
Poisson  parameter  is  the  same  for  all  J  aerosol  samples  obtained  during  the  sampling  time 
T  =  Jt.  The  primary  goal  is  to  detect  the  pre.sence  of  the  bio-target  signal  in  a  given 
realization  of  the  received  count  data  {F{j),  j  =  1,  2, . . . ,  J}  (i.e.,  in  the  J  aerosol  samples 
collected  during  the  sampling  time  T). 


Since  the  bio-target  signal  may  or  may  not  be  present  in  the  received  data  sequence 
{ F(j)}^Lii  this  leads  to  the  following  binary  hypothesis  testing  problem  involving  the  choice 
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between  the  two  hypotheses: 

^i:A«(j)  =  S(i)  +  AW.  (2) 


According  to  the  classical  theory  of  detection  [4],  it  is  desirable  to  implement  a  threshold 
signal  detection  receiver  which  realizes  the  likelihood  ratio  test.  The  likelihood  ratio  test  is 
preferred  because  it  possesses  some  nice  asymptotic  (i.e.,  large-sample  size)  properties  such 
as  unbiasedness,  consistency,  and  constant  faJse-alarm  rate.  With  reference  to  Eq.  (2),  the 
likelihood  ratio  test  decides  Hi  (i.e.,  bio-target  present)  if 


p{{FU)}\Ho) 


>  Ao, 


where  Aq  is  some  threshold  value  and  p{{F(j)}\Hi)  is  the  likelihood  function  for  obtaining 
the  count  data  {F(j),j  =  1,2,...,J}  under  hypothesis  H  (i  =  0,1).  Hence,  to  detect 
the  bio-target  signal,  we  have  to  compare  the  detection  statistic  A  with  a  threshold  Ao, 
the  “decision  level”.  The  exceedance  of  this  decision  level  by  A  is  an  indication  of  the 
presence  of  the  bio-target  signal.  For  Poisson-distributed  count  data,  it  is  simpler  to  use 
the  log-likelihood  ratio  test  which  decides  Hi  if 

£  =  In  A  =  ln(p({F(j)}|Fi))  -  ln(p({F(j)}|Fo))  >  T, 
where  T  =  In  Aq  is  the  threshold  for  the  transformed  decision  rule. 


Because  the  count  data  {F{j),j  —  1,2,  are  a.s.sumed  to  be  independently 

distributed,  the  computation  of  the  log-likelihood  ratio  is  straightforward.  This  results  in 
the  following  test: 


In  A 

.7 


=E 

.7=1 


F(j)ln  14- 


m') 

A(0)  ) 


-s{j) 


>T. 


(3) 


The  threshold,  T,  is  selected  objectively  using  the  Neyman- Pearson  criterion  [4].  In  the 
Neyman-Pearson  criterion,  the  probability  of  false  alarm  in  constrained  to  an  acceptable 
value  and  under  this  constraint,  the  decision  rule  is  designed  to  maximize  the  probability  of 
detection.  The  false-alarm  probability,  PpA,  (level  of  the  test)  and  the  detection  probability, 
Pd,  (power  of  the  test)  are  defined  as 


Ffa  =  Pr{£  >  T\Ho}, 


(4a) 
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and 

Po  =  Pr  {jC>T\Hi},  (46) 

where  Pr{-}  denotes  the  “probability  that”. 

To  select  a  threshold  level  that  corresponds  to  some  prespecified  false-alarm  proba¬ 
bility  requires  a  knowledge  of  the  probabilistic  characteristics  of  the  detection  statistic  C. 
The  determination  of  T  using  the  exact  statistics  of  C  (which  is  described  by  a  multivariate 
Poisson  probability  function)  is  a  formidable  problem.  However,  when  the  number  of  counts 
in  a  sample  is  sufficiently  largely  (e.g.,  >  50),  the  Poisson  probability  function  for  the 

count  data  can  be  approximated  very  well  with  the  Gaussian  distribution.  In  view  of  this, 
we  apply  the  Lindelberg  condition  [5]  of  the  central  limit  theorem  to  the  sum  of  random 
variables  to  argue  that  the  log-likelihood  ratio,  C,  exhibited  in  Eq.  (3)  is  asymptotically 
Gaussian-distributed  under  both  Hq  and  Hi  in  accordance  to  the  following  description; 

C  ~  N{fj,Oy<^o)  under  JTo; 

£  ~  iV(/ii,(Ti)  under  iJi-  (5) 

Here,  iV(^,  ct)  denotes  the  Gaussian  dLstribution  with  mean  //  and  standard  deviation  cr . 
The  mean  and  variance  of  the  log-likelihood  ratio  under  Ho  and  Hi  can  be  shown  to  have 
the  following  forms: 

/io  =  E{C\Ho)  =  ^  ln(l  -b  -  Sij)]  ,  (6) 

^i  =  E(£|H-i)  =  E(£jJJo)  +  +  w)’ 

i=i 

s  Var(£|fl-o)  =  E  +  w) ' 

i=l 

and  ^ 

al  =  Var(£|H'i)  =  Var(£|H-o)  +  ^  S{j)  In^ll  +  ^).  (9) 

i=i 

Here,  E(-)  denotes  mathematical  expectation,  and  Var(-)  denotes  variance.  An  examination 
of  Eqs  (6)-(9)  shows  that  the  variance  of  C  under  Hi  approaches  the  variance  of  C  under 
Ho  as  J  oo  and  S{j)  ->■  0  (for  all  j)  at  a  rate  that  is  faster  than  the  convergence  of  the 
means  under  the  two  hypotheses.  Consequently,  in  this  limit,  the  decision  test  for  £  reduces 
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to  a  test  between  two  Gaussian  variates  with  dijBFerent  means  and  the  same  variance.  In 
this  case,  the  detection  performance  is  completely  embodied  in  the  detectability  parameter 

E(£|Hi)  -E(£[go) 

y'Var(£|iTo) 


Eqs  (5)-(9)  provide  the  basis  for  an  easily-evaluated  computation  for  PpA  and  Pd- 
It  is  now  a  straightforward  exercise  to  show  that  the  value  of  the  threshold,  T,  required  to 
maintain  a  desired  false-alarm  probability  of  a  (significance  level)  can  be  determined  from 
the  solution  of  the  following  equation  [cf.  Eqs  (4a)  and  (5)]: 


PpA  =  a  = 


T  -  //o' 


O' 


where 


$(a;)  = 


1 


exp 


(4) 


dt 


is  the  cumulative  standard  normal  distribution  function.  Hence, 


(lOo) 


(106) 


T  =  /^  -I-  (7o$  “  “)  • 


(11) 


For  the  specified  threshold,  T,  the  probability  of  detection,  Pd,  [cf.  Eqs  (4b)  and  (5)]  is 
calculated  as 

=  (12) 

It  is  interesting  to  note  that  for  the  stressing  case  where  the  signal  intensity  level  is  much 
less  than  the  ambient  background  clutter  intensity  level,  the  probability  of  detection,  Pd,  for 
a  fixed  false-alarm  probability,  PpA,  is  a  monotonically  increasing  function  of  the  quantity 


(13) 


which  we  take  to  be  the  definition  of  the  signal-to-noise  ratio  for  our  problem. 
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III.  DISTRIBUTED  SENSOR  DETECTION  AND  DATA  FUSION 


In  this  section,  the  theory  of  distributed  detection  of  BW  agents  is  developed.  To 
this  end,  let  us  consider  the  following  two-stage  hypothesis  testing/detection  procedure. 
In  the  first  stage,  the  i-th  local  detector  (i  =  1, 2, . . . ,  n)  measures  the  local  aerosol  count 
data  {F^{j),j  =  1, 2, . . . ,  Jj}  and  applies  the  likelihood  ratio  test  to  produce  a  local  binary 
detection  decision  Dj  (z  =  1, 2, . . . ,  n);  viz., 

^  _  r  0,  if  detector  i  decides  Hq  (bio-target  absent  locally) ; 

*  ~  ^  1,  if  detector  i  decides  Hi  (bio-target  present  locally). 

In  the  second  stage,  these  binary  detection  decisions  Di  {i  =  1,2, ...  ,n)  are  then  transmitted 
to  the  data  fusion  center  where  they  are  processed  through  an  appropriate  fusion  rule 
(combining  strategy)  to  produce  a  global  detection  decision.  Da,  for  the  presence  or  absence 
of  the  bio-target  (i.e.,  a  global  inference);  viz.. 


_  J  0,  decide  bio-target  absent,  i.e.,  global  decision  Hq', 

^  ~  (  1,  decide  bio- target  pre.sent,  i.e.,  global  decision  Hi. 

It  is  assumed  that  the  local  detection  decisions  are  mutually  statistically  independent,  i.e., 
there  is  no  direct  communication  between  the  local  detector  elements  of  the  distributed 
system — the  local  detectors  are  only  permitted  to  send  the  local  detection  decisions  directly 
to  the  data  fusion  center  (central  processor,  or  equivalently,  the  command  and  control 
center).  To  summarize,  the  fusion  center  yields  the  global  decision.  Da,  based  on  the 
received  decision  vector  containing  the  local  decisions,  viz.,  D  =  (Di,D2, . . . ,  D„).  The 
global  decision.  Da,  is  assumed  to  depend  only  on  the  local  decision  vector,  D,  and  not  on 
the  observations,  {F^{j),j  =  1,2,...,  J,},  at  the  individual  detectors. 


The  objective  now  is  to  determine  a  set  of  appropriate  fusion  or  decision  rules  at  the 
fusion  center  for  the  specification  of  a  global  decision  that  maximizes  the  system  detection 
probability  under  a  false-alarm  probability  constraint.  Hence,  the  goal  is  to  formulate 
the  distributed  Neyman- Pearson  detection  problem  and  then  to  derive  the  decision  rules. 
Towards  that  goal,  we  seek  to  assign  local  detection  thresholds  %  {i  =  l,2,...,n),  for 
the  local  (peripheral)  likelihood  ratio  detectors  in  the  system  so  that  the  global  detection 
probability,  Pd,  at  the  fusion  center  (power  of  the  .system)  is  maximized  while  keeping 
the  global  false-alarm  probability,  Pfa,  at  a  desired  level,  say  a  (level  of  the  test).  This 
constrained  maximization  problem  can  be  solved  by  application  of  Lagrange’s  multiplier 
method  which  reduces  in  the  present  case  to  the  maximization  of  the  following  Lagrange 
function  L: 


L(Ti,T2,  . . . ,%) 


=  Pd  (Ti,  T2,  ...,%) +X  [PpA{ri,r2, . 


(14) 
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where  A  is  an  undetermined  Lagrange  multiplier.  The  dependence  of  the  global  probabilities 
of  detection  and  false  alarm  on  the  local  detection  thresholds  has  been  explicitly 

indicated  in  Eq.  (14).  Hence,  the  decision  rules  that  maximize  Pq  for  a  given  value  of  I^a 
are  likelihood  ratio  tests  with  thresholds  to  be  determined. 

To  utilize  Eq.  (14)  to  determine  the  optimal  detector  structure  requires  explicit 
expressions  for  the  global  detection  probability,  Pq,  and  the  global  false-alarm  probability, 
PpA-  The  precise  form  of  these  global  probabilities  depends  on  the  fusion  rule  used  to 
combine  the  detection  decisions,  Du  at  the  data  fusion  center.  More  precisely,  the  fusion 
rule  is  completely  specified  by  the  following  transition  probabilities: 


P(Di,D2,. 

..,Dn) 

=  Pr{DG  =  l|A,P2,. 

•  • ,  D„}, 

(15a) 

G[Di,D2,  ■ 

=  Pr{DG  =  0|Di,D2,. 

(156) 

In  what  follows,  let  PpAi  and  Pq.  denote  the  probabilities  of  false  alarm  and  detection, 
respectively,  for  the  i-th  local  detector.  Obviously,  these  two  probabilities  are  determined 
by  the  choice  of  the  local  detection  thrasholds  %  [cf.  Eqs  (10)-(12)].  The  constraint  Pfa  =  « 
induces  constraints  on  the  individual  Pfaj  and  determines  the  level  at  each  of  the  local 
detectors.  Because  Pq  is  an  increasing  function  of  Pp>.  (i  =  to  obtain  the 

maximum  value  of  fb,  it  follows  that  each  local  detector  should  operate  at  its  maximum 
achievable  Pd^,  viz.,  use  the  most  powerful  test,  for  its  given  level  PpAr 

Given  the  form  of  the  transition  probabilities  defined  in  Eq.  (15),  expressions  for 
the  global  probabilities  can  be  determined  in  terms  of  the  transition  probabilities  and  the 
local  probabilities  of  detection  and  false  alarm  as  follows:  (1)  the  global  probability  of  false 
alarm,  PpA, 

Ppa  =  Pv{Dg  =  1\Ho]  =  Pr{DG  =  l\D}-Pv{D\Ho} 

Di^D2f—iDn 

n 

=  Pv{DG  =  l\DuD2,...,Dn]-'[[Pv{Di\Ho} 

Di,D2,—,D„  »=1 

=  Y,  F{DuD2,...,D„)-1[{1-Ppa,)-IIPfa/,  (16) 

■^0 

(2)  the  global  probability  of  detection,  Pd, 

Pd  =  Pr{DG  =  =  Y 

D\,D2,...,D„ 
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=  PT{DG  =  l\DuD2,...,Dn}-l[Pr{Di\Hi] 

Di,D2,...,D„  t=l 

=  ^  F{DuD2,...,Dn)-Y[{l-PDj)-llPD,-,  (17) 

Di,D2,—  tDn  Mo  M] 

and,  (3)  the  global  probability  of  miss,  Pm, 

Pm  =  1-Pd=  Y.  Pr{PG  =  0|D}.Pr{D|i?i} 

Di,D2,...,D„ 

Di,D2,--M^n  Mo  Ml 

In  Elqs  (16)-(18),  the  index  set  Mq  =  [i:  Di  —  O);  similarly,  the  index  set  Mi  =  {i  :  Di  = 
l}.  In  other  words,  Mq  is  the  set  of  all  local  BW  agent  detectors  that  decide  “bio-target 
absent”,  whereas  Mi  is  the  set  of  all  detectors  that  decide  “bio-target  present”. 

With  these  explicit  expressions  for  Pq  and  PpA,  we  are  now  in  a  position  to  maximize 
the  Lagrange  function  L  [cf.  Eq.  (14)]  with  respect  to  the  local  detection  thresholds  % 
(i  =  1, 2, . . . ,  n).  By  the  usual  rules  of  calculus,  this  is  accomplished  by  equating  the  partial 
derivative  of  L  [cf.  Eq.  (14)]  with  re.spect  to  %  to  zero  for  i  =  1, 2, . . . ,  n.  This  leads  to  a 
system  of  n  nonlinear  equations  in  (n-f  1)  unknowns,  namely,  the  local  detection  thresholds 
{%]  and  the  Lagrange  multiplier  A.  Another  independent  equation  results  by  equating  the 
partial  derivative  of  L  with  respect  to  A  equal  to  zero.  This,  as  usual,  simply  encodes  the 
constraint  equation,  i.e.,  Pfa  =  oc. 

IV.  SOME  FUSION  RULES  AND  PERFORMANCE  ANALYSIS 

It  is  not  possible  to  obtain  explicit  solutions  for  the  optimum  detection  thresholds 
%  without  a  knowledge  of  the  combining  strategy  at  the  fusion  center.  Hence,  the  local 
detection  thresholds  %  can  be  determined  using  Lagrange’s  multiplier  method  once  the 
fusion  rule  for  combining  the  local  detection  decisions  Di  {i  =  1,  2, . . . ,  n)  has  been  specified. 
As  an  illustration  of  the  technique,  we  consider  a  combining  strategy  that  is  based  on  the 
fc-out-of-n  logical  function;  viz.,  if  k  or  more  local  detectors  decide  Hi,  then  the  global 
detection  decision  is  to  decide  Hi.  As  a  simple  illustration  of  the  technique,  we  consider  a 
distributed  system  consisting  of  three  local  detectors.  For  this  simple  case,  we  will  provide 
explicit  results  for  the  detection  performance  of  the  distributed  system  based  on  the  k- 
out-of-n  fusion  rule.  For  a  distributed  system  of  three  detectors  (i.e.,  n  =  3),  we  need 
to  consider  the  l-out-of-3  rule  (OR  rule),  the  2-out-of-3  rule  (MAJORITY  rule),  and  the 
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3-out-of-3  rule  (AND  rule).  In  what  follows,  denote  the  mean  and 

standard  deviation  of  the  log- likelihood  ratio  statistic  measured  at  the  i-th  local  detector 
under  hypotheses  Ho  and  Hi,  respectively  [cf.  Eqs  (6)-(9)]. 

Case  1:  (OR  rule)  For  this  case,  the  components  of  the  binary- valued  state  vector 
D  =  (Di,  D2,  D3)  are  combined  at  the  data  fusion  center  using  the  Boolean  logic  operator 
OR.  Consequently,  the  associated  transition  probabilities  are  given  by 

F(D).Pr{D.  =  l|D)  =  {;. 

Or,  equivalently, 

G(5).Pr{fl.=0|D}  =  {j. 

Here,  |D  |  =  D,  is  the  Hamming  norm  of  the  binary- valued  vector  D.  In  other  words, 
the  data  fusion  center  decides  Hq  (bio-target  absent)  only  if  all  the  local  detectors  decide 
Ho. 


For  the  case  of  the  OR  rule,  the  Neyman-Pearson  decision  rule  is  most  easily  ex¬ 
pressed  as  the  minimization  of  the  global  probability  of  miss,  Pm,  [which  is  equivalent  to 
the  maximization  of  the  global  probability  of  detection,  Pd  =  (1  -  Pm)]  subject  to  the 
constraint  that  the  global  probability  of  false  alarm,  PpA,  be  maintained  at  a  desired  level  a 
(i.e.,  PpA  =  a).  In  view  of  Eqs  (10),  (12),  (14),  (15),  (16),  and  (18),  the  choice  of  the  local 
detection  thresholds,  %,  to  achieve  optimum  detection  performance  can  be  formulated  as 
the  minimization  of  the  following  Lagrange  function: 


L(ri,  Ta,  T3)  =  Pm  +  A  [(1  -  Ppa)  -  (l  -  a) 


=n« 


i=l 


(O' 


-bA 


-  (1  -  a) 


Taking  the  partial  derivatives  81/8%  {i  =  1,2,3)  and  8L/8\  and  equating  these  results 
to  zero,  leads  to  the  following  set  of  design  relationships  for  the  specification  of  the  local 
detection  thresholds  required  to  achieve  optimum  global  detection  performance: 


,(0 


+  A 


^(0 


(O' 


i=i 


% _ [H 


(O' 


(2  =  1,2,3),  (19a) 
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where  $'(2;)  =  d^{x)ldx  and 


=  (!-“)• 


(196) 


Case  2:  (MAJORITY  rule)  In  accordance  with  this  combining  or  voting  rule,  the 
collective  decision  produced  at  the  data  fusion  center  concerning  the  presence  or  absence  of 
the  bio-target  is  that  made  by  more  than  half  of  the  local  detectors  of  the  system  (i.e.,  by 
at  least  two  detectors  in  the  case  of  a  system  of  three  detectors).  For  this  decision  strategy, 
the  transition  probabilities  reduce  to 

F(i5).Pr{D„  =  l|i3}  =  {j. 


Now,  with  reference  to  Eq.  (14),  we  consider  the  maximization  of  the  following  Lagrange 
function: 


T(ri,  T2,  T3)  =  Pd  +  A  [i^FA  -  a] 

i=l  j=l  i=l 


+  A 


^(1  -  PfaJ  ■  n  S  -  a 


t=l 


3=1 

jY* 


i=l 


Taking  the  partial  derivative  of  L  with  respect  %  {i  =  1,2,3),  setting  the  result  to  zero, 
and  inserting  the  explicit  expressions  for  PpA^  and  Pd-  [cf.  Eqs  (10a)  and  (12)]  result  in  the 
following  system  of  equations  (i  =  1,2,3): 


■  E 

j,k=i 


3 

■  E 

j,k=l 


+  A 


1_$ 


^0) 


$ 


(20a) 
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This  equation,  when  supplemented  with  the  constraint, 


E* 

1=1 


(i) 


a, 


n  1-* 


j=i 


^0) 


or 


+  E  1-^ 


1=1 


%zj3 

ji) 


or 


=  a,  (206) 


permits  the  calculation  of  the  local  detection  thresholds  for  optimum  global  detection  per¬ 
formance. 

Case  3:  (AND  rule)  For  this  combining  rule,  the  data  fusion  center  decides  Hi  only 
if  all  the  local  detectors  decide  Hi ,  so  that  the  transition  probabilities  are  given  by 

F(D)  =Pr{DG--l\D}  =  ll' 

^  ^  ^  \  0,  otherwise; 


and 


G(S).Pr{Z).=0|D}={j. 

Substituting  these  transition  probabilities  into  Elq.  (14),  the  Lagrange  function  reduces  to 

L(Ti, 72,73)  =  Fb  +  a[^7Va  -  a 
=nPb,+A 


i=l 


Lt=l 


Now  equating  dLjdFi  (i  =  1,2,3)  to  zero  leads  to 


0  =  -- 


(O' 


n  i-$ 


%  Fj 
^0) 


(.O' 


.;=i 


+  A 


(O' 


.(0 


.(0 


n 


j=i 


..(0 


(O' 


(i  =  1,2,3).  (21a) 


'We  need  to  solve  these  equations  subject  to  the  constraint 

,(0\\ 


=  a. 


(216) 


Since  an  analytical  solution  of  the  system  of  nonlinear  algebraic  equations  for  the 
optimal  local  detection  thresholds  is  not  possible,  it  is  necessary  to  use  a  numerical  method. 
Although  the  system  of  four  equations  in  the  four  unknowns  %  (i  =  1,2,3)  and  A  is  highly 
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nonlinear,  this  system  of  equations  can,  nevertheless,  be  readily  solved  using  a  quasi-Newton 
method  that  utilizes  an  appropriate  line-search  algorithm,  as  described  by  Dennis  and 
Schnabel  [6]. 

In  order  to  provide  some  concrete  results,  we  proceed  to  compare  the  achievable 
global  detection  performance  of  a  distributed  system  of  three  local  detectors  with  that 
provided  by  a  single  detector.  In  so  doing,  we  will  demonstrate  the  improvement  of  the 
detection  performance  of  the  distributed  .system  of  detectors  over  that  provided  by  a  single 
detector.  For  the  numerical  results,  we  consider  a  constant  target  aerosol  signal  S{j)  = 
So  for  j  =  1, 2, . . . ,  J,  where  J  is  the  number  of  aerosol  samples  used  in  the  detection 
procedure.  For  the  simulations  presented  here,  a  constant  bio-target  signal  So  was  adjusted 
to  provide  the  desired  signal-to-noise  ratio  defined  in  accordance  to  Eq.  (13).  We  have  used 
a  background  aerosol  clutter  count  with  an  intensity  =  10.  The  detection  performance 
for  the  distributed  .system  using  the  three  fusion  rules  (i.e.,  OR,  MAJORITY,  and  AND 
rules)  and  for  a  single  detector  are  compared  in  terms  of  the  probability  of  detection  (Pd) 
determined  as  a  function  of  the  signal-to-noi.se  ratio  (SNR)  with  both  the  sample  size,  T,  (or, 
equivalently,  detection  time)  and  the  probability  of  false  alarm,  Pfa,  as  implicit  parameters. 
The  numerical  results  are  presented  in  Figs  1-9  for  different  values  of  these  parameters. 
Each  of  these  figures  contains  four  detection  performance  curves,  corresponding  to  either 
the  single  detector  or  the  distributed  system  for  one  of  the  three  combining  strategies. 
These  curves  have  been  calculated  for  observation  sample  sizes  of  J  =  10,  20,  and,  30  and 
comsidered  in  combination  with  probabilities  of  false  alarm  of  Pfa  =  10~^,  10“^,  and  10 

The  following  conclusions  can  be  drawn  fi-om  Figs  1-9.  Firstly,  since  the  detec¬ 
tion  performance  curve  for  a  single  detector  is  always  the  lowest  curve  in  these  figures,  the 
improved  performance  achieved  with  the  distributed  system  is  clearly  noticeable.  We  see 
that  for  a  given  SNR,  sample  size,  and  false-alarm  probability,  the  distributed  system  always 
performs  better  (in  the  sense  of  a  higher  detection  probability)  than  the  single  detector.  Sec¬ 
ondly,  as  to  be  expected,  the  detection  performance  of  the  single  and  distributed  detection 
system  monotonically  increases  as  the  observation  sample  size,  J ,  and/ or  the  false-alarm 
probability,  Pfa,  increases.  Thirdly,  the  detection  performances  of  the  distributed  system 
for  the  three  fusion  rules  considered  are  comparable  for  the  larger  observation  sample  sizes 
(e.g.,  for  J  =  30)  and/or  for  the  larger  false-alarm  probabilities  (e.g.,  for  Pfa  =  10“^). 
Fourthly,  for  the  smaller  sample  sizes  and/or  false-alarm  probabilities,  the  detection  perfor¬ 
mance  of  the  distributed  system  is  more  strongly  dependent  on  the  fusion  rule  (combining 
strategy)  used  to  combine  the  local  detection  decisions  at  the  fusion  center.  Even  so,  it  is 
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not  possible  to  rank  the  fusion  rules  with  respect  to  detection  performance  over  the  entire 
range  of  SNRs  considered,  since  it  is  apparent  that  no  single  fusion  rule  is  optimal  over 
the  complete  SNR  regime.  However,  a  perusal  of  Figs  1-9  shows  that  the  MAJORITY 
rule  generally  provides  a  better  detection  performance  for  a  given  set  of  conditions  (e.g., 
fixed  SNR,  sample  size,  and  false-alarm  probability)  than  either  the  OR  or  AND  rule  for 
SNRs  ranging  from  about  —25  dB  to  -5  dB.  Furthermore,  the  detection  performance  of 
the  AND  rule  is  superior  to  that  of  either  the  MAJORITY  or  the  OR  rule  in  the  lower 
SNR  regime  (e.g.,  less  than  about  —30  dB).  Finally,  the  crossover  point  (in  dB)  at  which 
the  performance  of  one  fusion  rule  is  superior  to  that  of  another  (in  the  sense  of  better 
detection  performance)  is  a  complex  function  of  both  the  observation  sample  size  and  the 
false-alarm  probability. 

A  clearer  picture  of  the  relative  detection  performances  of  the  single  detector  and 
the  distributed  system  for  the  three  fusion  rules  can  be  obtained  if  we  define  a  minimum 
detectable  signal  (MDS)  in  terms  of  the  SNR  (in  dB)  required  to  yield  a  detection  probability 
Pd  =  0.5  for  a  fixed  false-alaiTn  probability,  Pfa-  The  MDS  has  been  computed  as  a 
function  of  Pfa  values  from  10'"*  to  10“^.  These  results  are  displayed  in  Fig.  10  for  the 
single  detector  and  for  a  distributed  system  using  one  of  the  three  fusion  rules  considered 
above.  An  observation  sample  size  of  J  =  20  has  been  used  to  evaluate  the  MDS.  It  is  clear 
from  Fig.  10  that  the  MDS  can  be  used  to  rank  the  detection  performance  for  the  single 
detector  and  for  the  distributed  system  using  one  of  the  three  fusion  rules.  Observe  that  the 
minimum  detection  performance  of  the  distributed  system  is  superior  to  that  of  the  single 
detector.  Furthermore,  we  note  that  in  the  case  of  a  distributed  system,  the  MAJORITY 
rule  clearly  exhibits  superior  minimum  detection  perfonnance  relative  to  either  the  AND  or 
OR  rule.  Indeed,  over  the  range  of  Pfa  values  considered,  we  determine  a  value  of  about  3.5 
dB  improvement  in  minimum  signal  detectability  for  a  distributed  system  employing  the 
MAJORITY  rule  over  the  single  detector.  The  MAJORITY  rule  for  the  distributed  system 
exhibits  an  MDS  performance  improvement  of  w  0.5  dB  and  w  3.0  dB  over  that  provided 
by  the  AND  rule  and  OR  rule,  respectively. 

V.  PLACEMENT  OF  DETECTORS  IN  A  DISTRIBUTED  SYSTEM 


The  concentration  field  arising  from  a  BW  agent  source  is  inhomogeneous,  non¬ 
isotropic,  and  time  dependent  as  a  result  of  the  source  configuration  and  the  non-stationarity 
and  inhomogeneity  of  the  underlying  atmospheric  (environmental)  flow.  Because  the  BW 
agent  cloud  concentration  is  neither  homogeneous  nor  isotropic,  it  is  important  to  consider 
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the  proper  placement  of  the  detectors  in  a  distributed  system  in  order  to  maximize  detection 
performance. 

As  the  BW  agent  puflp-cloud  passes  over  the  target  area,  it  is  necessary  to  ensure  that 
the  individual  detectors  in  the  distributed  system  be  located  so  that  the  detectors  have  a 
high  probability  of  being  inside  the  “core”  of  the  BW  agent  cloud  where  the  concentrations 
(and  dosages)  are  high.  Clearly,  the  proper  choice  of  the  number  of  samples,  J ,  to  use  in  the 
detection  process,  and  the  choice  of  thresholds  %  required  to  optimize  the  detection  perfor¬ 
mance  Pd  of  the  system  for  a  given  SNR  [cf.  Eq.  (13)],  while  ensuring  simultaneously  that 
a  certain  target  probability  of  false  alarm  Pp/^  is  never  exceeded,  is  an  extremely  complex 
problem.  The  choice  of  the  appropriate  detection  and  false-alarm  probabilities  in  a  par¬ 
ticular  operational  context  is,  of  course,  a  question  for  military  commanders  to  decide  and 
one  on  which  physical  science  has  no  bearing.  But  it  needs  to  be  stressed  that  (in  general) 
operational  standards  cannot  be  based  on  Pd  =  1  with  Pp\  =  0;  it  can  never  be  absolutely 
guaranteed  that  the  distributed  system  will  provide  perfect  detection  performance  under 
any  particular  operating  .scenario.  Once  an  operating  point  is  decided  upon  (i.e.,  once  Pd 
and  Pfa  have  been  chosen),  the  choice  of  J  is  set  by  the  maximum  allowable  mean  time 
to  detection  (viz.,  when  a  bio- target  in  some  BW  agent  is  encountered  in  the  target  area, 
the  objective  is  to  minimize  the  mean  time  to  detection  of  the  “threatening”  bio-target  so 
that  appropriate  countermeasures  can  be  taken).  The  maximum  allowable  mean  time  to 
detection  can  be  specified  a  priori  based  on  operational  considerations  and  this  value  can 
be  used  subsequently  to  determine  an  upper  bound  for  J.  With  this  choice  of  J ,  and  given 
an  operating  point  specified  by  Pd  and  Pfa,  receiver  operating  characteristic  (ROC)  curves 
similar  to  those  exhibited  in  Figs  1-9  can  be  used  to  determine  the  SNR  required  to  achieve 
this  operating  point  for  the  distributed  system.  In  some  situations,  such  as  long-range 
stand-off  BW  agent  releases  where  the  agent  concentration  at  the  detectors  is  small,  the 
SNR  required  to  achieve  a  particular  operating  point  in  the  distributed  system  may  not 
exist  and  the  current  methodology  will  fail  to  provide  a  solution  in  this  case.  This  simply 
implies  that  the  desired  system  performance  cannot  be  achieved  in  this  case,  but  even  in 
this  case,  the  methodology  yields  useful  information  in  the  sense  that  it  provides  the  user 
with  a  specific  prediction  on  when  the  detection  system  will  fail. 

The  selection  of  J  and  SNR  now  permits  the  determination  of  the  minimum  dosage 
that  needs  to  be  acquired  by  the  individual  BW  agent  detectors  in  order  achieve  the  target 
detection  performance.  In  particular,  the  minimum  required  dosage  can  be  determined 
from  Eq.  (13).  Naturally,  the  minimum  dosage  level  required  for  detection  depends  on  the 
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background  aerosol  clutter  count  whose  intensity  is  embodied  in  the  parameter  It 

is  assumed  that  the  background  has  been  characterized  a  priori  (viz.,  the  noise  intensity 
parameter  for  the  operating  environment  is  assumed  to  be  known).  The  minimum 
dosage  at  each  BW  agent  detector  required  to  achieve  the  target  detection  performance 
can  be  determined  as  follows  from  Eq.  (13): 

.7 

=  J]S(j)  =  JAo-SNR,  (22) 

i=i 

where  the  dosage  Dmin  must  be  achieved  over  the  total  sampling  time  T  =  Jr.  To  ensure 
that  the  individual  BW  agent  detectors  receive  at  least  Pmin  during  the  passage  of  a  BW 
agent  cloud  over  the  target  area,  the  detectors  must  be  located  in  a  region  that  will  be 
covered  with  the  threshold  dosage  "Dmin  or  higher  with  a  high  probability. 

The  problem  of  the  proper  placement  of  BW  agent  detectors  in  order  to  achieve  a 
certain  global  probability  of  detection  for  the  distributed  system  now  reduces  to  the  problem 
of  the  determination  of  the  locus  of  points  in  the  target  region  where  the  dosage  is  equal  to 
or  greater  than  "Dmini  with  some  specified  probability.  In  other  words,  we  need  to  deterndne 
the  statistical  properties  of  the  dosage  and  the  fluctuating  cloud  concentration  from  which 
it  is  derived.  To  that  purpose,  we  need  to  investigate  the  statistical  properties  of  the  partial 
dosage  accumulated  over  a  sampling  (detection)  time  interval  T  =  Jr,  on  exposure  of  a 
detector  located  at  x  =  (x,  y,  z)  to  the  instantaneous  concentration  Xs(x,  t)  of  the  dispersing 
BW  agent  cloud:  ^ 

'D{x,T)=[  Xs(x,t)dt.  (23) 

Jo 

In  Eq.  (23),  it  has  been  assumed  without  any  fundamental  loss  of  generality  that  the  lower 
limit  t  =  0  corresponds  to  the  time  at  which  the  detection  interval  began.  The  time  origin 
t  =  0  is  then  assumed  to  be  reset  for  each  new  detection  interval  of  duration  T.  Because 
Xs{x.,t)  is  a  random  function,  V{x,T)  mast  be  a  random  (stochastic)  integral.  Hence,  for 
a  fixed  point  x  and  a  fixed  sampling  time  T,  P(x,  T)  is  a  random  variable. 

Naturally,  the  statistical  properties  of  V  must  be  determined  by  the  probabilistic 
characteristics  of  the  underlying  cloud  concentration  fluctuations.  The  general  probabilistic 
behavior  of  the  random  variable  V  can  be  expressed  by  its  probability  density  function 
(PDF);  and  the  essential  statistical  quantities  that  give  information  about  the  shape  of  this 
function  are  the  values  of  the  integral  moments.  Hence,  if  we  define  a  PDF  for  the  partial 
dosage  V  at  position  x  over  a  sampling  time  T  as 

/(D;x,r)dD  =  Pr{D<P(x,r)  <D  +  dD},  (24) 
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the  fc-th  order  integral  moment  of  V  is  given  by 

AOO 

fXk{V-T)  =  {V’^iT))=  D*/(D;x,T)dD 

Jo 

dti  dt2  ...  dtk,  fc  €  N;  (25a) 

or, 

poo 

fj.[iV-,T)  =  {iV{T)-fii)^)=  /  iD-tiyffiD-,^,T)dD 

Jo 

{x's(ti)x's{t2) . . .  x'sitk))  dh  dh...  dtk,  fc  e  N,  (256) 

where  the  primes  on  a  fluctuating  quantity  denote  departures  from  the  mean  [e.g.,  x'si^)  = 
{Xs{t)  -  (Xs(i)))])  3.nd  (...)  denotes  an  emsemble  average  (or,  mathematical  expectation). 
Here,  the  dosage  is  implicitly  assumed  to  be  determined  at  a  fixed  point  so  the  explicit 
x-dependence  in  Eqs  (25a)  and  (25b)  has  been  suppressed. 

Note  that  the  fc-th  order  integral  moments  of  the  dosage  V  require  complete  informa¬ 
tion  on  the  joint  fc-time  point  moments  of  the  fluctuating  cloud  concentration  Xs{t)-  Hence, 
in  order  to  obtain  the  statistical  properties  of  T>,  we  require  a  complete  knowledge  of  all 
the  higher-order  multi-time  statistics  of  the  fluctuating  cloud  concentration  field.  Unfortu¬ 
nately,  very  little  information  is  currently  available  concerning  the  higher-order,  multi-time 
point  statistics  of  a  cloud  concentration  field.  In  particular,  only  one-  and  two-time  point 
statistics  of  the  cloud  concentration  appear  to  have  been  measured  or  modelled.  This  infor¬ 
mation  permits  the  determination  of  the  two  lowest-order  moments  of  the  dosage;  namely, 
the  mean  and  variance  of  V  which  we  will  consider  briefly  below. 

In  view  of  Equation  (25a)  evaluated  for  fc  =  1,  the  ensemble-mean  dosage  for  a 
statistically  non-stationary  cloud  (puff)  release  reduces  to  the  following  form: 

T 

fz,{V;T)  =  {V{T))  =  {xs{t))dt.  (26) 

Equation  (26)  implies  that  the  ensemble-mean  dasage  can  be  determined  from  only  a  knowl¬ 
edge  of  the  ensemble- mean  concentration  of  the  BW  agent  cloud. 

However,  because  the  dosage  is  random  owing  to  the  fluctuating  cloud  concentration, 
it  is  necessary  (at  the  very  least)  to  also  estimate  the  variance  of  the  dosage  in  order  to 
provide  the  simplest  measure  of  the  variability  that  can  be  expected  in  this  quantity  between 
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two  realizations  of  the  dispersion.  Towards  this  objective,  evaluation  of  Equation  (25b)  for 
k  =  2  leads  to  the  desired  expression  for  the  dosage  variance; 


^  lo  L 

/  PiXs-,ti,s)ds;  (27a) 


where 


p{Xs;t,r) 


{x's('t)x'A'^)) 

{Xs^(i)) 


(276) 


is  the  autocorrelation  function  for  the  Xs{t)  random  process;  t  and  r  are  diffusion  times. 
Hence,  (t  —  r)  corresponds  to  the  correlation  lag  time.  We  note  that  the  concentration  of  a 
diffusing  cloud  measured  at  a  fixed  point  in  space  does  not  constitute  a  stationary  process 
with  the  result  that  its  autocorrelation  function  is  a  function  of  two  diffusion  times,  t  and 
r  [or,  equivalently,  of  a  diffusion  time  t  and  a  lag  time  (t  —  r)]. 


We  now  express  the  dosage  variance  as  follows: 

rT  .  _  ,  fti 


(V''‘(T))  =2  j  /  p{x.;h,‘)ds 

where  T;^(t)  is  the  integral  time  scale  for  the  Xs{t)  process  at  the  diffusion  time  t: 

Tx(.t)  =  f  P(Xf;  s)  ds. 

Jo 


(28) 


(29) 


Thus,  in  order  to  calculate  the  dosage  variance,  we  need  an  estimate  of  the  integral  time 
scale  of  concentration  fluctuations.  A  simple  estimate  of  T^{t)  has  been  provided  by  Sykes 
[7]  who  used  Gifford’s  [8]  meandering  plume  model  to  provide  a  simple  scheme  for  esti¬ 
mating  of  plume  concentration  from  a  knowledge  of  the  Eulerian  velocity  time  scale. 
This  analysis  showed  that  the  effect  of  intermittency  is  to  reduce  the  integral  time  scale  of 
concentration  fluctuations  from  the  Eulerian  wind  fluctuation  time  scale,  and  also  to  intro¬ 
duce  a  logarithmic  correction  factor  that  depends  on  the  relative  fluctuation  intensity  (i.e., 
standard  deviation  of  concentration  divided  by  mean  concentration)  of  the  concentration 
fluctuations. 


Given  a  model  for  the  ensemble-mean  concentration,  concentration  variance,  and 
integral  time  scale  of  concentration  fluctuations,  it  is  now  po.ssible  to  determine  the  mean 
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and  variance  of  the  dosage  using  Eqs  (26)  and  (28),  respectively.  Unfortunately,  the  dosage 
V{T)  cannot  really  be  adequately  described  in  terms  of  the  ensemble-mean  dosage  and 
dosage  variance  without  also  specifying  the  form  of  the  dosage  PDF.  Our  current  knowl¬ 
edge  of  cloud  concentration  statistics  do  not  enable  us  to  predict  an  exact  form  for  the 
dosage  PDF  (i.e.,  the  latter  would  effectively  require  explicit  knowledge  of  the  multi- time 
point,  higher-order  statistics  of  the  cloud  concentration  field).  In  consequence,  to  make 
further  progress  the  only  recourse  we  have  available  at  present  is  to  compare  dosage  data, 
derived  from  measured  atmospheric  concentration  time  .series,  to  some  common  functional 
forms  for  model  PDFs.  To  that  end,  it  should  be  mentioned  that  analysis  of  a  large  number 
of  plume  concentration  time  series  measured  in  the  cooperative  Concentration  Fluctuation 
Experiments  (CONFLUX)  project  [9-11]  has  shown  that  the  dosage  PDF  can  be  approxi¬ 
mated  reasonably  well  by  the  clipped-normal  PDF  form  for  sampling  times  T  greater  than 
about  5  s  [12].  In  particular,  when  the  clipped-normal  PDF  is  fitted  to  the  dosage  data 
derived  from  the  fluctuating  plume  concentration,  it  appears  to  represent  most  features  of 
the  statistics  of  the  data  (e.g.,  various  higher-order  moments  of  the  dosage,  quantiles  of  the 
dosage,  etc.)  quite  well,  at  least  for  exposure  (averaging)  times  T  >  5  s.  In  this  regard,  it 
should  perhaps  be  noted  that  the  effect  of  smoothing  (i.e.,  by  increasing  the  averaging  or 
sampling  time  over  which  the  dosage  is  determined)  generally  tends  to  improve  the  fit  of 
the  clipped- normal  distribution  to  the  random  dosage  compared  to  that  of  the  underlying 
fluctuating  concentration  process  from  which  it  was  derived. 

In  consequence,  the  dosage  PDF  will  be  assumed  to  be  adequately  modelled  using 
a  clipped-normal  PDF  which  is  defined  as  follows: 


where 


^1  -f-  erf 


J^\ 
\f2ao  ) 


(31) 


and  ^lD  and  ao  are  location  and  scale  parameters,  respectively.  In  Eq.  (31),  erf(-)  denotes 
the  error  function.  In  the  clipped- normal  PDF,  a  Gaus.sian  PDF  is  u.sed  to  model  the 
continuous  part  of  the  distribution  of  D  for  P  >  0,  with  the  probability  of  negative  dosages 
being  transferred  into  a  delta  function  at  =  0  to  represent  the  intermittency  effect.  In 
the  clipped-normal  PDF,  the  intermittency  factor  7  is  uniquely  determined  as  a  function  of 
//£)/cr£).  Also,  the  two  parameters  fip  and  ap,  which  completely  define  the  clipped-normal 
PDF,  can  be  uniquely  determined  from  the  mean  dosage,  (D(T)),  and  dosage  variance. 
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(X>'^(T)),  using  the  following  two  relationships: 

and 

{V'\T))  =  -{V{T)f  +  GDi+tJ^D{V{T)),  (33) 

where  7  is  defined  in  Equation  (31). 

We  may  now  exploit  the  information  contained  in  the  statistical  description  of  the 
partial  dosage  for  the  proper  placement  of  individual  BW  agent  detectors  in  a  distributed 
system  in  order  to  satisfy  certain  constraints  on  the  global  detection  probability  for  a  fixed 
false-alarm  rate  with  a  prescribed  allowable  maximum  mean  time  to  detection.  These 
constraints  can  be  met  if  the  individual  BW  agent  detectors  are  placed  at  points  in  the 
target  region  where  the  dosage  is  greater  than  or  equal  to  a  “critical”  or  “threshold”  dosage 
prescribed  in  Eq.  (22).  In  other  words,  for  a  prescribed  “critical”  or  “threshold”  dosage 
Pmin,  we  are  interested  in  the  determination  of  the  points  in  the  target  region  where  the 
dosage  received  in  a  particular  sampling  (detection)  interval  at  least  equal  to  X>min,  with 
some  specified  probability.  The  goal  is  to  place  the  individual  BW  agent  detectors  of  the 
distributed  system  within  an  area  of  the  target  region  that  will  be  well  covered  to  at  least 
the  critical  dosage  ©min  over  a  given  sampling  (detection)  time. 

In  view  of  Eq.  (30),  the  value  of  F{D\-k,T)  where 

poo 

F{D-,k,T)=  /(D';x,r)dD',  (34) 

Jd 

is  the  probability  that  any  one  realization  of  the  ensemble  of  the  partial  dosages,  accu¬ 
mulated  over  sampling  time  T  at  x,  is  greater  than  D.  Alternatively,  it  is  the  proportion 
of  realizations  in  the  whole  ensemble  of  BW  agent  clouds  for  which  the  partial  dosage  V 
received  in  a  particular  sampling  (detection)  time  T  and  at  location  x  is  greater  than  D. 
The  equation 

F{D-^,T)  =  a,  (35) 

where  a  is  a  fixed  con.stant  between  zero  and  one  [i.e.,  a  €  (0, 1)]  is,  for  a  particular  sampling 
time  T,  a  surface  at  every  point  of  which  there  is  a  probability  a  that  the  partial  dosage  © 
obtained  in  duration  T  exceeds  D.  For  every  point  outside  the  surface,  the  probability  that 
the  partial  dosage  accumulated  in  a  sampling  time  T  exceeds  D  is  less  than  a.  Furthermore, 
for  fixed  T,  each  point  in  space  is  associated  with  just  one  value  of  a  so  that  it  lies  on  a 
unique  surface  of  the  family  of  surfaces  that  is  specified  by  Eq.  (35).  Each  of  these  surfaces 
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intersects  the  plane  =  k  {k  is  a  constant)  in  a  curve  or  contour.  Here,  Zg  refers  to  the 
height  of  the  detectors  above  the  ground  surface.  To  ensure  a  high  probability  of  detection 
for  the  distributed  system,  the  individual  BW  agent  detectors  should  be  placed  within  that 
contour  for  which  the  dosage  level  is  at  least  the  “critical”  or  “threshold”  dosage  X>min 
prescribed  in  Eq.  (22) . 

VI.  SUMMARY  AND  CONCLUSIONS 


In  this  report,  we  have  presented  a  methodology  for  the  design  and  placement  of 
BW  agent  detectors  in  a  distributed  system.  The  first  step  of  the  methodology  involves 
the  modelling  of  the  detection  performance  of  the  individual  detectors.  In  the  regard,  we 
have  focussed  on  the  generalized  likelihood  ratio  algorithm  based  on  a  Poisson  probability 
model  for  binary  hypothesis  testing.  This  model  is  used  to  describe  the  counting  statistics 
associated  with  the  detection  of  a  bio-target  signal  embedded  in  a  common  aerosol  clutter 
background.  By  using  asymptotic  normal  theory,  the  approximate  distribution  of  the  de¬ 
tection  statistic  (viz.,  the  log-likelihood  ratio  in  this  case)  has  been  derived.  We  showed 
how  this  can  be  used  to  predict  the  probability  of  detection  for  a  prescribed  probability  of 
false  alarm  as  a  function  of  the  threshold  level,  the  signal-to- noise  ratio  of  the  detection 
environment,  and  the  total  number  of  samples  used  for  the  detection. 

The  second  step  of  the  methodology  concerns  the  design  and  characterization  of  the 
combining  strategy  that  will  be  used  to  provide  the  final  global  detection  decision.  For 
a  given  BW  agent  remote  detection  scenario,  the  detection  and  false-alarm  probabilities 
for  each  local  detector  in  the  distributed  system  can  be  expressed  explicitly  in  terms  of 
its  thresholds  and  operating  SNRs.  These  expressions  can  be  substituted  into  a  set  of 
n  equations  for  the  optimum  thresholds  derived  in  Section  III.  This  results  in  a  set  of  n 
equations  in  the  n  unknown  thresholds.  These  equations  involve  a  common  undetermined 
(Lagrange)  multiplier  A.  Together  with  the  overall  target  false-alarm  probability  constraint 
Pfa  =  where  Pfa  is  given  in  Eq.  (16),  the  resulting  design  equations  may  be  solved 
numerically  for  the  optimum  thresholds  for  a  given  combining  strategy  at  the  data  fusion 
center.  The  overall  system  detection  performance,  Pd,  can  then  be  evaluated.  When  this 
procedure  has  been  carried  out  for  all  the  combining  strategies  at  the  data  fusion  center 
and  for  a  given  set  of  operating  SNRs,  we  then  choose  that  strategy  which  results  in  the 
best  performance.  The  constraint  imposed  by  the  operating  SNRs  can  be  achieved  by  the 
correct  placement  of  the  local  detectors,  which  constitutes  the  third  and  final  step  of  the 
methodology. 
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Explicit  results  are  provided  for  the  specific  case  of  a  distributed  system  of  three 
detectors  with  the  fusion  rule  at  the  data  fusion  center  based  on  the  fc-out-of-n  logical 
function.  For  the  case  of  a  distributed  system  of  three  detectors  (i.e.,  n  =  3),  this  logical 
function  provides  the  OR  rule  {k  =  1),  the  MAJORITY  rule  {k  =  2),  and  the  AND  rule 
{k  =  3).  Numerical  results  of  the  global  detection  performance  of  the  distributed  system  for 
these  three  fusion  rules  are  provided  and  this  performance  is  compared  to  that  achievable 
with  a  single  detector.  The  numerical  simulations  and  results  indicate  that,  for  a  given  set 
of  operating  conditions  (e.g.,  signal-to-noise  ratio,  sample  size,  and  false-alarm  probability), 
the  distributed  system  always  exhibits  better  detection  performance  (e.g.,  in  the  sense  of 
a  higher  detection  probability)  than  that  afforded  by  a  single  detector.  Furthermore,  we 
showed  that  for  the  distributed  system,  no  single  fusion  rule  appears  to  provide  an  optimal 
detection  performance  over  the  entire  range  of  SNRs.  It  has  been  observed  that  the  order 
of  detection  performance  of  the  distributed  detector  .system,  under  the  various  fusion  rules 
used,  is  a  complex  function  of  the  signal-to-noise  ratio,  the  number  of  observations  used  in 
the  detection  process,  and  the  probability  of  false  alarm.  With  an  increase  in  the  sample 
size,  J,  (or,  equivalently,  an  increase  in  the  sampling  time)  and/or  an  increase  in  the 
tolerable  false-alarm  probability,  the  detection  performance  of  both  the  single  detector  and 
the  distributed  system  improved.  Furthermore,  with  an  increase  in  the  sample  size  and/or 
allowed  false-alarm  probability,  the  differences  in  performance  of  the  various  combining 
strategies  (fusion  rules)  narrow. 

The  third  step  of  the  methodology  involves  the  correct  placement  of  the  individual 
detectors  in  the  distributed  system  over  the  target  area.  The  selection  of  sample  size,  J, 
which  in  turn  is  determined  by  the  maximum  allowable  mean  time  to  detection,  and  of 
the  operating  signal-to-noise  ratio  permits  the  determination  of  the  minimum  dosage  level 
that  needs  to  be  acquired  by  the  individual  (peripheral)  BW  agent  detectors  in  order  that 
the  overall  system  achieve  a  target  global  detection  performance.  In  particular,  the  target 
detection  performance  for  the  distributed  system  can  be  met  if  the  individual  BW  agent 
detectors  can  be  placed  at  points  in  the  target  region  where  the  dosage  level  reaches  at  least 
a  “critical”  or  “threshold”  dosage  Pmin  [cf-  Eq.  (22)].  This  requires  the  determination  of 
the  points  in  the  target  region  where  the  dosage  received  in  a  particular  sampling  (detec¬ 
tion)  time  is  equal  to  or  higher  than  ‘Dminj  with  some  specified  probability.  In  consequence, 
a  model  for  the  ensemble  mean  dosage  and  dosage  variance  has  been  derived  that  allows 
these  quantities  to  be  determined  from  an  explicit  knowledge  of  the  mean  concentration, 
concentration  variance,  and  an  integral  time  scale  of  concentration  fluctuations  for  the  dis¬ 
persing  BW  agent  cloud.  It  is  important  to  note  that  there  currently  exists  a  practical 
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probabilistic  model  for  turbulent  diffusion  [i.e.,  the  Second-Order  Closure  Integrated  PUFF 
(SCIPUFF)  model  [13]  developed  by  Aeronautics  Research  Associates  Princeton  (ARAP) 
under  the  sponsorship  of  Defence  Special  Weapons  Agency  (DSWA)]  that  is  capable  of  pro¬ 
viding  estimates  for  the  mean  concentration,  concentration  variance,  and  an  integral  time 
scale  of  concentration  fluctuations.  The  incorporation  of  the  proposed  methodology  for  de¬ 
tector  placement  within  SCIPUFF  will  permit  then  a  completely  probabilistic  or  statistical 
approach  to  the  correct  placement  of  individual  BW  agent  detectors  in  a  distributed  system 
for  a  given  BW  agent  release  scenario. 
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Figure  1 


Detector  performance  comparison  (expressed  as  the  probability  of  detection,  fb,  versus 
SNR)  of  a  distributed  system  composed  of  three  detectors  using  the  OR,  MAJORITY, 
and  AND  fusion  rules  at  the  central  processor  and  of  a  single  detector  for  J  =  10  and 
PpA  =  10“^. 
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Figure  2 

Detector  performance  comparison  (expressed  as  the  probability  of  detection,  Pd,  versus 
SNR)  of  a  distributed  system  composed  of  three  detectors  using  the  OR,  MAJORITY, 
and  AND  fusion  rules  at  the  central  processor  and  of  a  single  detector  for  J  =  10  and 
PPA  =  10-3. 
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Figure  3 

Detector  performance  comparison  (expressed  as  the  probability  of  detection,  Fb,  versus 
SNR)  of  a  distributed  system  composed  of  three  detectors  using  the  OR,  MAJORITY, 
and  AND  fusion  rules  at  the  central  processor  and  of  a  single  detector  for  J  =  10  and 
PpA  =  10~^. 
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Figure  4 

Detector  performance  comparison  (expressed  as  the  probability  of  detection,  Pdi  versus 
SNR)  of  a  distributed  system  composed  of  three  detectors  using  the  OR,  MAJORITY, 
and  AND  fusion  rules  at  the  central  processor  and  of  a  single  detector  for  J  =  20  and 

PpA  =  10“^. 
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Figure  5 


Detector  performance  comparison  (expressed  as  the  probability  of  detection,  Pq,  versus 
SNR)  of  a  distributed  system  composed  of  three  detectors  using  the  OR,  MAJORITY, 
and  AND  fusion  rules  at  the  central  processor  and  of  a  single  detector  for  J  =  20  and 
Pfa  =  10-3. 
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Figure  6 

Detector  performance  comparison  (expressed  as  the  probability  of  detection,  fb,  versus 
SNR)  of  a  distributed  system  composed  of  three  detectors  using  the  OR,  MAJORITY, 
and  AND  fusion  rules  at  the  central  processor  and  of  a  single  detector  for  J  =  20  and 
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Figure  7 


Detector  performance  comparison  (expressed  as  the  probability  of  detection,  Pdi  versus 
SNR)  of  a  distributed  system  composed  of  three  detectors  using  the  OR,  MAJORITY, 
and  AND  fusion  rules  at  the  central  processor  and  of  a  single  detector  for  J  =  30  and 
Pfa  =  10“^. 
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Figure  8 

Detector  performance  comparison  (expressed  as  the  probability  of  detection,  Pq,  versus 
SNR)  of  a  distributed  system  composed  of  three  detectors  using  the  OR,  MAJORITY, 
and  AND  fusion  rules  at  the  central  processor  and  of  a  single  detector  for  J  —  30  and 
PfA  =  10-3. 
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Figure  9 

Detector  performance  comparison  (expressed  as  the  probability  of  detection,  Pd,  versus 
SNR)  of  a  distributed  system  composed  of  three  detectors  using  the  OR,  MAJORITY, 
and  AND  fusion  rules  at  the  central  processor  and  of  a  single  detector  for  J  =  30  and 
PpA  = 
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Figure  10 

The  minimum  detectable  signal  (MDS)  versus  the  probability  of  false  alarm,  Pfa>  for  ^ 
distributed  system  composed  of  three  detectors  using  the  OR,  MAJORITY,  and  AND 
fusion  rules  at  the  central  processor  and  of  a  single  detector. 


DRES-SR-669 


UNCLASSIFIED 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  FORM 
{highest  classificBtion  of  Title,  Abstract,  Keywords) 


DOCUMENT  CONTROL  DATA 

(Sacurity  elatsification  of  titia,  body  of  abatraet  and  Indaxing  annotation  must  ba  antarad  whan  tha  ovarall  doeumant  it  clattifiad) 

1.  ORIGINATOR  (the  name  and  address  of  the  organization  preparing  the  document 

2.  SECURITY  CLASSIFICATION 

Organizations  for  whom  the  document  was  prepared,  e.g.  Establishmem  sponsoring 

(overall  security  classification  of  the  document 

a  contractor's  report  or  tasking  agency,  are  entered  in  section  8.) 

including  special  warning  terms  if  applicable) 

Defence  Research  Establishment  Suffield 

P.O.  Box  4000 

UNCLASSIFIED 

Medicine  Hat,  Alberta  TIA  8K6 

3.  TITLE  (the  complete  document  title  as  indicated  on  the  title  page.  Its  classification  should  be  indicated  by  the  appropriate 
abbreviation  {S,C,R  or  U)  in  parentheses  after  the  title.) 


Methodology  for  the  Design  and  Optimal  Placement  of  Point  Detectors  in  a  Distributed  Detection  System 
for  Remote  Defence  Against  Biological  Warfare  Agent  Releases  (U) 


4.  AUTHORS  (Last  name,  first  name,  middle  initial.  If  military,  show 

Yee,  Eugene  C. 

rank,  e.g.  Ooe,  Maj.  John  E.) 

5.  DATE  OF  PUBLICATION  (month  and  year  of  publication  of 

6a.  NO.  OF  PAGES  (total 

6b.  NO.  OF  REFS  (total  cited  in 

document) 

containing  information.  Include 

document) 

February,  1998 

Annexes,  Appendices,  etc.) 

13 

35 

6.  DESCRIPTIVE  NOTES  (the  category  of  the  document,  e.g.  technical  report,  technical  note  or  memorandum.  If  appropriate,  enter  the  type  of 
report,  e.g.  interim,  progress,  summary,  annual  or  final.  Give  the  inclusive  dates  when  a  specific  reporting  period  is  covered.) 


Technical  Report  (Final)  [1  January,  1997  to  31  December,  1997] 


8.  SPONSORING  ACTIVITY  (the  name  of  the  department  project  office  or  laboratory  sponsoring  the  research  and  development  Include  the 
address.) 


Defence  Research  Establishment  Suffield 
P.O.  Box  4000,  Medicine  Hat,  Alberta  TIA  8K6 


9a.  PROJECT  OR  GRANT  NO.  (if  appropriate,  the  applicable  research 
and  development  project  or  grant  number  under  which  the  document 
was  written.  Please  specify  whether  project  or  grant) 

PCN  No.  051SP 

9b.  CONTRACT  NO.  (if  appropriate,  the  applicable  number  under 
which  the  document  was  written) 

10a.  ORIGINATOR'S  DOCUMENT  NUMBER  (the  official  document 
number  by  which  the  document  is  identified  by  the  originating 
activity.  This  number  must  be  unique  to  this  document) 

SR-669 

10b.  OTHER  DOCUMENT  NOS.  (Any  other  numbers  which  may 
be  assigned  this  document  either  by  the  originator  or  by  the 
sponsor) 

1  1.  DOCUMENT  AVAILABILITY  (any  limitations  on  further  dissemination  of  the  document,  other  than  those  imposed  by  security  classification) 

(X)  Unlimited  distribution 

(  )  Distribution  limited  to  defence  departments  and  defence  contractors;  further  distribution  only  as  approved 

(  )  Distribution  limited  to  defence  departments  and  Canadian  defence  contractors;  further  distribution  only  as  approved 

(  )  Distribution  limited  to  government  departments  and  agencies;  further  distribution  only  as  approved 

(  )  Distribution  limited  to  defence  departments;  further  distribution  only  as  approved 

(  )  Other  (please  specify): 

1 2.  DOCUMENT  ANNOUNCEMENT  (any  limitation  to  the  bibliographic  announcement  of  this  document  This  will  normally  correspond  to 

the  Document  Availabilty  (11).  However,  where  further  distribution  (beyond  the  audience  specified  in  11)  is  possible,  a  wider 
announcement  audience  may  be  selected.) 

Unlimited 

UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  FORM 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  FORM 

13.  ABSTRACT  (  a  brief  and  factual  summary  of  the  document  It  may  also  appear  elsewhere  in  the  body  of  the  document  itself.  It  is  highly 
desirable  that  the  abstract  of  classified  documents  be  imclassified.  Each  paragraph  of  the  abstract  shall  begin  with  an  indication  of  the 
security  classification  of  the  information  in  the  paragraph  (unless  die  documem  itself  is  unclassified!  represented  as  (S).  (0.  (Rl,  or  (U). 

It  is  not  necessary  to  include  here  abstracts  in  both  offical  languages  unless  the  text  is  bilingual). 

This  report  deals  with  the  design  and  performance  evaluation  of  a  distributed  detection  system  for  a 
dispersing  biological  warfare  (BW)  cloud  embedded  in  the  natural  aerosol  component  of  the  background 
air.  The  distributed  system  employs  a  number  of  physically  separated  BW  agent  point  detectors  (sentries) 
located  within  some  target  region  and  a  data  fusion  center  that  provides  the  final  decision  as  to  the  presence 
or  absence  of  the  bio-target  by  combining  the  individual  localized  decisions  from  the  various  point  detectors 
using  a  prespecified  combining  strategy.  In  this  system,  each  detector  implements  a  generalized  likelihood 
ratio  test  on  its  own  localized  observations  to  test  for  the  presence  or  absence  of  a  bio-target.  These  localized 
detection  decisions  are  then  transmitted  to  a  data  fusion  center  where  they  are  logically  combined  to  yield  a 
global  detection  decision  for  the  distributed  system.  The  optimization  of  the  global  detection  performance 
of  the  distributed  system  is  derived  by  application  of  the  Lagrange  multipliers  method,  whereby  the  global 
probability  of  detection  for  the  system  is  maximized  subject  to  the  constraint  that  the  global  probability  of 
false  alarm  is  maintained  at  a  prespecified  level  (viz.,  at  a  constant  and  tolerable  false-alarm  rate).  A  number 
of  combining  strategies  are  investigated  in  order  to  determine  some  overall  system  optimality  for  detection.  It 
is  found  that  the  optimal  detection  strategy  for  each  individual  BW  agent  detector  in  the  distributed  system 
depends  on  the  strategy  of  all  the  other  detectors,  as  well  as  on  the  structure  of  the  data  fusion  center. 
Explicit  numerical  results,  in  the  form  of  the  probability  of  detection  versus  signal-to-noise  ratio  for  several 
preassigned  false-alarm  probabilities  and  sample  sizes  [which  determine  the  maximum  allowable  mean  time 
to  detection],  are  presented  for  the  case  of  a  distributed  system  consisting  of  three  BW  agent  detectors.  The 
results  are  compared  to  the  detection  performance  that  is  achievable  using  only  a  single  detector.  Finally,  the 
problem  of  the  correct  placement  of  the  individual  BW  agent  detectors  in  the  distributed  system  in  order  to 
achieve  a  specified  detection  performance  is  addressed.  In  particular,  a  method  for  the  determination  of  the 
dosage  probabilities  at  fixed  points  in  the  target  region  is  developed  and  used  to  place  BW  agent  detectors 
in  a  distributed  system. 

14.  KEYWORDS,  DESCRIPTORS  or  IDENTIFIERS  (technically  meaningful  terms  or  short  phrases  that  characterize  a  document  and  could  be 
helpful  in  cataloguing  the  document  They  should  be  selected  so  that  no  security  classification  is  required.  Idemifiers,  such  as  equipment 
model  designation,  trade  name,  military  project  code  name,  geographic  location  may  also  be  included.  If  possible  keywords  should  be  selected 
from  a  published  thesaurus,  e.g.  Thesaurus  of  Engineering  and  Scientific  Terms  (TEST)  and  that  thesaurus-idemtfied.  If  it  is  not  possible  to 
select  Indexing  terms  which  are  Unclassified,  the  classification  of  each  should  be  indicated  as  with  the  title.) 


Detection  Algorithm  and  Performance 
Probability  of  Detection  and  False- Alarm 
Signal  Processing  and  Analysis 
Bio-aerosol  Detection 
Placement  of  Detectors 
Distributed  Detection 
Sensor/Data  Fusion 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  FORM 


